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ABSTRACT 

A set of cloud retrieval algorithms developed for CERES and applied to MODIS data have been adapted to analyze 
other satellite imager data in near-real time. The cloud products, including single-layer cloud amount, top and base 
height, optical depth, phase, effective particle size, and liquid and ice water paths, are being retrieved from GOES- 
10/11/12, MTSAT-1R, FY-2C, and Meteosat imager data as well as from MODIS. A comprehensive system to 
normalize the calibrations to MODIS has been implemented to maximize consistency in the products across platforms. 
Estimates of surface and top-of-atmosphere broadband radiative fluxes are also provided. Multilayered cloud properties 
are retrieved from GOES- 12, Meteosat, and MODIS data. Native pixel resolution analyses are performed over selected 
domains, while reduced sampling is used for full-disk retrievals. Tools have been developed for matching the pixel- 
level results with instrumented surface sites and active sensor satellites. The calibrations, methods, examples of the 
products, and comparisons with the ICESat GLAS lidar are discussed. These products are currently being used for 
aircraft icing diagnoses, numerical weather modeling assimilation, and atmospheric radiation research and have 
potential for use in many other applications. 
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1. INTRODUCTION 

With the availability of well-calibrated polar-orbiting imagers such as the Moderate Resolution Imaging 
Spectroradiometer (MODIS) 1 on the Terra and Aqua satellites, faster computers, new multispectral imagers on 
geostationary satellites (GEOsats), and spaceborne lidars and radars on the Cloud- Aerosol Lidar and Infrared Pathfinder 
Satellite Observations (CALIPSO) 2 , the Ice, Cloud, and Elevation Satellite (ICESat) 3 , and CloudSat 4 , it is now possible 
to provide, in near-real time, more accurate cloud and radiation products than ever before. There remain, however, 
many obstacles to effectively utilize those resources to actually produce cloud and radiation properties for use in 
nowcasting and numerical weather prediction models over the globe. This paper describes a prototype system developed 
to analyze global geostationary satellite data at many scales for both research and potential operational applications. 

2. DATA 


2.1 Satellite data 

Imager radiance data from the Geostationary Operational Environmental Satellites (GOES-E and GOES-W), and the 
European Meteosat Second Generation, the Chinese Fengyun (FY-2C), and the Japanese Multi-Functional Transport 
Satellite (MTSAT-1R) satellites functionally provide coverage at all latitudes equatorward of ~60° latitude. Terra and 
Aqua or the operational polar-orbiting meteorological satellites can provide coverage at higher latitudes. While all of 
these satellites have imagers with five or more channels, not all are useful for cloud monitoring. As seen in Table 1, all 
of the referenced satellites have three similar spectral channels that are useful for cloud detection and retrievals, and six 
of the imagers hold a fourth relevant channel (split-window, 12.0 pm) in common. Four of the imagers have a CCL- 
slicing channel that is also used for cloud retrievals. The nominal resolutions common to all channels on a given imager 
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Table 1. Satellite imager channels used in common for cloud detection and retrievals. 


Satellite 

Sub-satellite 
Longitude or 
Equator Crossing 
Time 

Common 

(visible) 

Central or radiance equivalent channel wavelength, pm 


resolution, 

km 

Visible 

(VIS) 

Shortwave-IR 

(SIR) 

Infrared 

(IR) 

Split-window 

(SW) 

co 2 

GOES- 11 

135°W 

4(1) 

0.651 

3.91 

10.7 

12.0 

N/A 

GOES-12 

75°W 

4(1) 

0.651 

3.92 

10.7 

N/A 

13.3 

Meteosat 

0° 

3(1) 

0.640 

3.92 

10.8 

12.0 

13.4 

FY-2C 

105°E 

5(1.25) 

0.724 

3.79 

10.8 

12.0 

N/A 

MTSAT- 1R 

140°E 

5(1.25) 

0.724 

3.79 

10.8 

12.0 

N/A 

Aqua 

1330 LT 

1 (0.25) 

0.646 

3.79 

10.8 

12.0 

13.3 

Terra 

1030 LT 

1 (0.25) 

0.646 

3.78 

10.8 

12.0 

13.3 


vary from 1 to 5 km. Higher resolution visible (VIS) channels are available on all of the imagers. Full-disk data from the 
GEOSats are available from, at least, hourly ( Meteosat , MTSAT-R, FY-2C) to 3-hourly (GOES). The Meteosat full- 
disk data are available at higher temporal resolution, while nearly full-disk images can be constructed hourly from the 
GOES subsection images. Higher temporal resolution images of subsections are available from all of the GEOSats. The 
polar-orbiting satellites sample a given area 6-8 times per day in the polar latitudes. In this paper, the MODIS data are 
only used in reference to calibration; their use in near-real time global cloud analyses is considered in future studies. 

2.2 Satellite calibrations 

Each GEOSat imager channel is calibrated against the corresponding Terra or Aqua MODIS channel using a variety of 
techniques 5 " 7 , unless the channel spectral differences are too large (e.g., 3.79 vs. 3.92 pm) or sampling is too limited. 
They are also cross-calibrated against each other 8 . If the differences are small, the nominal calibrations are used. The 
satellite calibrations are performed each month and the coefficients are updated online for the VIS channels. The Terra 
and Aqua VIS channels were found to be very stable and their measured reflectances differ by only 1% for a given 
scene 9 . Figure 1 shows the Meteosat - 9 and MTSAT-1R VIS channel calibrations against the Terra and Aqua MODIS 
VIS channels, respectively, for June 2007. The GEOSat data are given in brightness or squared brightness counts 
(CNT). The Meteosat-9 data (Fig. la) are tightly correlated with the Terra radiances, especially at the low end. The 
slopes of these linear regression fits vary little from month to month and exhibit a slight upward trend (Fig. lb), 
suggesting some degradation in the Meteosat-9 VIS response. Scatter between the MTSAT squared counts and Aqua 
radiances is much greater, even at the low end. Although, the MTSAT- 1R data are also available in 10-bit counts that 
are supposed to be linear in radiance, the scatterplots and fits are very similar. The spread of the MTSAT counts at the 
low end is a consistent, known problem that still seeks a solution. The FY -2C VIS channel yields similar results, but 
with more scatter. The correlations between the MODIS and GOES VIS channels are similar to those in Fig. la. 
Problems in accurately calibrating the MTSAT and FY-2C VIS channels introduce uncertainty in the retrievals. 
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Fig. 1. Correlations and regression fits between (a) Meteosat-9 and Terra MODIS VIS channels and (c) MTSAT-1R and Aqua 
MODIS VIS channels, June 2007. Standard linear fit in black, principal component fit in green, reversed axis linear fit in blue, and 
force fit to zero or known space count in red. (b) slopes of linear regression fits between Meteosat - 9 and Terra MODIS VIS channels. 
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Fig. 2. Same as Fig. 1, except for IR channels and (c) Term versus FY-2C. 

The IR channel calibrations show similar differences between the reference and GEOSats. The Meteosat- 9 IR 
brightness temperatures are highly correlated with and show little scatter relative to the Term values (Fig. 2a) but are, 
on average, 0.8 K greater than the Terra values during the day. At night (not shown), the difference is 0.6 K. This small 
day-night discrepancy is reflected in the small differences between the monthly slopes (Fig. 2b). The fits between the 
Terra and other GEOSat IR channels are similarly correlated, except for FY-2C. Both GOES IR channels have a 0.3 K 
day-night discrepancy relative to the Terra IR channel, but the average bias is close to zero. No day-night difference is 
evident in the MTS AT- Terra bias, which is less than 0.2 K. The relationships between the FY-2C and Terra IR 
temperatures vary from month to month. As seen Fig. 2c, the greatest discrepancies occur at lower temperatures. In this 
case, the Terra values are greater than their FY-2C counterparts for temperatures, T < 250 K. During other months, 
there is agreement at the low end, while, at other times, the low Terra temperatures are less than the FY-2C values. 
While part of this variation may be attributable to navigation errors in the FY-2C data, it is mirrored in cross- 
calibrations with MTS AT and Meteosat (not shown) indicating some significant stability uncertainties in the FY-2C IR 
data. The comparisons between the Terra split-window channels and those on each GEOSat are very much like those 
for the IR channels. 

The GEOSat SIR channels are also relatively well behaved, except for the FY-2C. The correlations among the GOES- 
11, GOES-12, and Meteosat- 9 are all very tight (Fig. 3b, c), although the GOES-12 measures higher temperatures than 
Meteosat- 9. The MTSAT SIR temperatures are highly correlated with their Aqua counterparts (Fig. 3a) and the mean 
day and night Aqua-WYSAJ SIR temperature differences are 0.7 and 0.4 K, respectively. Day-night differences in the 
GOES-11/12 and GOES -\2IMeteosat-9 biases are around +0.1 K. The FY-2C SIR channel is generally well behaved 
during the daylight hours, yielding good correlations and month-to-month consistency in the calibration fits against 
Aqua data, despite significant temperature differences. The average daytime slope and offsets arel.085 and -24.85 K. 



Fig. 3. Correlations and regression fits between daytime SIR data from (a) MTSAT-1R and Aqua MODIS, August 2006 and on (b) 
GOES-1 1 and GOES-12, and (c) GOES-12 and Meteosat-9, June 2007. Different fits explained as in Fig. 1. 
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Fig. 4. Nighttime FY-2C images, 1800 UTC, 13 August 2008 (a) 1R and (b) SIR-IR temperature difference, (c) Correlations and 


regression fits between nighttime SIR data from (a) FY-2C and MTSAT-1R. Different fits explained as in Fig. 1. 


At night, there are a number of erratic irregularities in the FY-2C SIR channel that include an east-west gradient and 
strange, varying striping patterns, as well as wholesale shifts in the gain. Figure 4 shows an example of two of the 
irregularities and the typical impact on intercalibration. The 1800 UTC, 13 August 2008 IR image (Fig. 4a) shows 
apparently normal behavior, while the corresponding SIR-IR brightness temperature difference (BTD) image (Fig. 4b) 
reveals a double cross near the middle and a large area in the east with values exceeding 20 K. Other large areas appear 
to be perfectly normal. This type of behavior is seen nearly every day between 1500 and 1800 UTC rendering those 
images nearly useless for inclusion in the calibrations and retrievals, but the behavior at other hours is more normal. 


Shortwave (SW) and longwave (LW) fluxes at the top of the atmosphere (TOA) are computed from the VIS and IR 
data, respectively, using narrow-to-broadband conversion functions 10 . The coefficients for those functions have been 
updated 11 using narrowband data matched to measured broadband fluxes from the Clouds and the Earth’s Radiant 
Energy System (CERES) project 12 . 


2.3 Auxiliary data 

Surface temperatures and profiles of temperature, humidity, and ozone are interpolated to a given location and time 
from the 6-hourly, 1.25° NOAA Global Forecast System (GFS) analyses and forecasts. Clear-sky albedos, surface 
emissivities, surface types, and snow and ice cover are estimated as for CERES 13 . 


3. METHODOLOGIES 


The pixel-level data are processed in a tile format. Each tile consists of all pixels in a 0.5° or 1° region, depending on 
the domain. Currently, the GEOSat data full-disk images are sampled to achieve a nominal resolution between 8 and 10 
km, depending on the satellite. The surface temperatures, clear-sky radiances, surface types, and spectral atmospheric 
corrections are computed once for each tile and then used to classify each pixel as clear or cloudy 13 . Cloud properties 
are retrieved for each cloudy pixel using the VIS-IR-SIR-SWC Technique (VISST) during the day and the SIR-IR-SWC 
Technique (SIST) at night 14 . Daytime is defined as all times when the solar zenith angle (SZA) is less than 82°. At 
night, the SIST is modified to use the CCE channel instead of the SWC for GOES-12. In that instance, the CO 2 and IR 
channels are used to assign the cloud height and optical depth, while the IR and SIR channels are used to estimate the 
particle size. At night, the cloud particle sizes, optical depths, and water paths are only estimated for optically thin 
clouds (optical depths less than 3 or so). For optically thick clouds, default values are used at night. 

Table 2 lists the retrieved pixel-level cloud properties and radiation parameters. The near-infrared (NIR; 1.62 or 2.13 
pm) reflectance and water vapor (WV; 6.7 pm) temperatures are included for the satellites having those channels. The 
scene ID or cloud mask indicates whether the pixel is cloudy, snow-free clear, or snow-covered clear. Aircraft icing 
potential is derived from the liquid water path, effective temperature, and effective droplet size 13 . The viewing zenith 
(VZA) and relative azimuth (RAA) angles are included along with the SZA for each pixel. 
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Table 2. Output parameters from GEOSat analyses. 


Pixel-Level Products - all pixels 

VIS reflectance 

N1R reflectance 

BTDfSIR - IR) 

WV temperature 

IR temperature 

BTD(IR-SWC) 

BTD(IR-C0 2 ) 

scene ID (cloud mask) 

TOA shortwave albedo 

TOA longwave flux 

icing potential 

SZA, VZA, RAA 

latitude 

longitude 




Pixel-Level Products 

i - cloudy pixels only 


phase 

optical depth 

droplet/crystal size 

liquid/ice water path 

effective temperature 

effective pressure 

top pressure 

bottom pressure 

IR emissivity 

effective height 

top height 

bottom height 

Gridded Products 

total cloud percentage 

total cloud property means 

total cloud property SDs 

mean SZA, VZA, RAA 

level cloud percentages 

level cloud property means 

mid cloud property means 

high cloud property means 

VIS reflectance 

IR temperature 

SIR temperature 

clear-sky IR temperature 

clear-sky VIS reflectance 

clear-area skin temperature 

mean column RH 

TOA shortwave albedo 

TOA longwave flux 

surface SW flux up 

surface SW flux down 

surface LW flux up 

surface LW flux down 





The pixel-level products are also averaged over the tiles to produce gridded properties. The cloud amount and the 
standard deviations (SD) are also provided for the total cloud properties. Similar averages are computed for the low (< 2 
km), middle 2-6 km) and high (> 6 km) cloud layers. Several other parameters are also computed for each tile 
including the surface radiative fluxes 11 , the surface skin temperature, VIS reflectance, and IR temperature for clear 
areas, and the mean relative humidity (RH) derived from the GFS data for the column above the surface. Currently, not 
all gridded products are available for all domains. 

The results of the analyses and the means for displaying and comparing them are described by Palikonda et al. 16 . These 
data and the tools can be accessed at http://www-angler.larc.nasa.gov/satimage/products.html. 

4. RESULTS 

Examples of the combined full-disk products are shown in Fig. 5 for 1800 UTC, 13 June 2008. Except for Meteosat-9, 
the analyses are confined to latitudes equatorward of 60° latitude. The pseudocolor RGB image (Fig. 5a) shows the 
terminator while providing some continuity between sunlit and dark portions of the globe. Areas of sunglint are evident, 
centered around 115°W, 80°W, and 30°W near 15°N. Except for GOES-11 and 12, the demarcation of the different 
satellites is easy to see, especially for FY-2C. The odd striping in the FY-2C SIR channel appears in this image also. 
Cloud-top height (Fig. 5b) also shows some satellite boundary lines, but for the Meteosat-9, the terminator boundary 
shows up in some areas (e.g., Africa) as a change in cloud height. Overall, the cloud-top heights are relatively consistent 
from one satellite to the next. The calibration problems for FY-2C seem to mostly affect the heights of the very thin 
high clouds. The optical depths (OD) in Fig. 5c and effective droplet radii (Re) in Fig. 5d also show the terminator lines 
because of the limitations of using thermal channels alone. The optical depths in the daylight sector range up to much 
larger values than over the dark regions. The effective droplet radii in the sunlit areas are very reasonable and show 
patterns similar to those seen from MODIS retrievals 14 , relatively small values near the coasts and much larger droplets 
in open marine areas. Currently, the full-disk results are only provided online in separate formats. A combined product 
such as that shown here will be available in the future when some of the satellite calibration and terminator problems 
have been resolved satisfactorily. 

The derived cloud and radiative properties have been compared against a variety of surface, aircraft, and other satellite 
datasets to assess their accuracies 11 ' 17 " 19 . Evaluation of the retrieved cloud properties is continuing as more reference 
data become available from various satellites and surface sites. The October 2003 cloud layering product from the 
ICESat Geoscience Laser Altimeter System 20 (GLAS) were matched with the GOES cloud products over the contiguous 
USA using two methods: one assigns GLAS footprints to individual GOES pixels and the second assigns the GLAS 
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Fig. 5. Combined imagery and cloud properties from GOES-1 1, GOES-12, Meteosat-9, FY-2C, and MTSAT-1R, 1800 UTC, 13 June 
2008. (a) RGB image created using combination of VIS, 1R, and SIR data. 
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Fig. 6. Daytime cloud amounts and differences over the contiguous USA domain, October 2003. (a) Combined mean 0.5° GOES- 
10/12 cloud amounts, (b) Differences between ICESat GLAS and matched GOES 0.5° cloud amounts. 

footprints to a given 0.5° region. Results for the second method are shown in Fig. 6 for daytime cloud fraction. Cloud 
amount differences of +5% predominate over areas with large and small mean cloud amounts. Comparisons of cloud 
amounts and heights yielded similar results for both approaches. During the daytime, the mean GOES cloud amounts 
were 52% for both methods compared to 55% from GLAS. At night, the GLAS average cloud amount was 64% because 
of its increased sensitivity to optically thin clouds, those with optical depths less than ~ 0.3. The GOES cloud amounts 
dropped to 51%. These results are similar other comparisons indicating the passive measurements detect few clouds 
with OD < 0.3. 

On average, the GOES cloud-top heights were 4.8 km during the day compared to 5.2 km from GLAS. At night, the 
bias increased from -0.4 to -1.9 km, primarily because the SIST was not using the CO 2 channel for GOES-12 at that 
time. Multi-layered clouds tend to cause an underestimate of cloud-top height. If only single-layer clouds are 
considered, the VISST tends to overestimate cloud-top height while the SIST still underestimates cloud height. A more 
comprehensive study 21 over a surface site in Oklahoma showed better agreement in single-layer cloud heights, perhaps 
because of more sampling and the use of GOES10/11 data instead of GOES-12. Initial comparisons of GLAS heights 
derived from GOES- 12 using a new IR-COj multi-layered cloud retrieval method 21 show significant improvement in the 
heights. 


5. CONCLUDING REMARKS 

A comprehensive system has been developed for producing nearly global cloud products in near-real time for a variety 
of applications including numerical weather model assimilation and nowcasting. The results from regional and full-disk 
analyses are now available online for use by the community. While many improvements can be made to the system, the 
greatest remaining obstacles to producing a continuously reliable product are consistent calibrations across all platforms 
and new techniques to address issues around the terminator. Retrievals over snow-covered surfaces are also potentially 
biased because of the decrease in contrast between the clouds and snow. Availability of 1.6-pm channels on Meteosat 
and the future GOES-R will help alleviate that problem. Integration of the new IR-CO 2 multilayered cloud retrieval 
method, currently applied to GOES- 12 and Meteosat data and applicable to GOES-R, into this system will improve the 
accuracy of the cloud-top heights. Validation studies using GLAS, CALIPSO, CloudSat and surface site data will 
continue and the results will be used to update the analysis codes. 
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